N87-18870 


ON COMPUTING EIGENSOLUTION SENSITIVITY DATA USING FREE 
VIBRATION SOLUTIONS 


B. P. Wang 

Department of Mechanical Engineering 
The University of Texas at Arlington 
Arlington, Texas 


SUMMARY 


A simplified method of computing eigensolution sensitivity derivatives in 
structural dynamics is developed in this paper. It is shown that if the elements 
of stiffness and mass matrices associated with a design variable are homogeneous 
functions of that design variable, then eigenvalue derivatives can be computed from 
element strain and kinetic energies. Furthermore, if cross-mode energies are known, 
eigensolution derivatives of modified systems can be computed approximately using 
assume mode reanalysis formulation. A ten bar truss example is used to illustrate 
the present formulations. 


INTRODUCTION 


The usefulness of eigensolution sensitivity derivatives in structural dynamics 
research is well known. The sensitivity data can be used for approximate 
reanalysis, analytical model improvement, assessment of design trend as well as 
structural optimization with eigenvalue constraints. When applied to larger discrete 
structural models, these applications typically require long and expensive computer 
runs and usually the predominate contributor to the computing time was the 
calculation of derivatives. Thus efficient eigensolution sensitivity analysis 
procedures would be very useful in structural dynamic research. It is the purpose 
of this paper to develop, under certain conditions, efficient eigensolution analysis 
procedures using free vibration data. 

The equations for computing derivatives of eigenvalues and eigenvectors for free 
vibration of undamped structures were known for a long time. Only recently have 
these methods been implemented in some general-purpose finite-element programs. In 
this paper, a simple method is developed which can be used to compute the eigenvalue 
derivatives for a large class of problems by exploiting the similarity between the 
equations for eigenvalue derivatives and element strain and kinetic energies. 
Furthermore, if the cross-mode element energy data are available, the approximate 
eigenvector derivatives can also be computed using a truncated modal expansion 
expression. The approximate second derivatives of eigenvalues can then be computed. 
Additionally, with the cross-mode strain energy data, the eigenvalue derivatives of a 
modified structure can be computed using assumed mode reanalysis formulation. 
Numerical examples will be presented to illustrate the various formulations. 
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EIGENSOLUTION SENSITIVITY IN STRUCTURAL DYNAMICS 


The general problem is to compute the rate of change (or derivatives) of 
eignevalues and eigenvectors with respect to design variables for the following 
generalized eigenvalue problem in structural dynamics. 

K 4> = A M 4 (1) 

Much research has addressed this problem in the past two decades. A 
comprehensive survey of literature can be found in a recent paper by Adelman and 
Haftka [1]. The equations for first order eigenvalue and eigenvector derivatives 
as well as second order eigenvalue derivatives are summarized below: 
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Second Derivative of Eigenvalues: 
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Note that in the above equations, the mode shapes are normalized to unit generalized 
mass, i.e. 


♦i M4> £ = 1.0 


234 


For the eigenvector derivatives, if less than full modes are used, Eq. (2) is an 
approximate expression. These will lead to approximate second order derivatives of 
eigenvalues . 

The above equations have been developed in the literature for some time. For 
example, Equations (2) and (3) can be found in Fox and Kapoor [2] and Eq. (7) was 
reported by Miura and Schmit [3]. It should be noted in passing that there are some 
recently developed algebraic methods [4-5] which can be used to compute eigenvector 
derivatives without using modal expansions. 

The difficulty of applying the aforementioned equations appears to be the cal- 
culation of derivatives of stiffness and mass matrices with respect to design vari- 
ables. In the next section it will be shown that under certain assumptions, we can 
circumvent the calculation of 3K/3x^ and 3M/3x£ in implementing these above equations. 

SIMPLIFYING ASSUMPTIONS 


In general, 


the system stiffness and mass matrices in Eq. (1) can be written as 


ND 

K = K + I K. (9) 
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where 

K = contribution to stiffness matrix due to structural elements that are to 

c 

remain constant during the design process. 

= contribution to mass matrix due to the masses of the unchanged elements as 
well as nonstructural masses. 


K^,M^ = contributions to stiffness and mass matrices respectively due to elements 
controlled by design variable x^. 

To develop simplified efficient methods for eigensolution sensitivity analysis, 

the elements of the matrices K. and M. are assumed to be homogeneous functions of 

design variables. That is thematrices K. and M. have the form 
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where K. and M. are constant matrices. Furthermore, define non-dimensional design 
i i 

parameters 
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Similarly, 


M. 
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where and M. n are stiffness and mass matrices due to design variable x. at its 


iO 


nominal value x. A . 

iO 

Based on the above assumptions, the derivatives of stiffness and mass matrices 
with respective to design variables can be computed readily: 
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Similarly, we can derive 
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It should be noted that at the nominal design, OL = OK = 1. With these simplifica- 
tions, the eigenvalue derivatives can be computed readily. 


RATE OF CHANGE OF EIGENVALUES 


Using (17) and (18) with o * a » 1, the eigenvalue derivative, eq. (2) becomes 
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Then Eq. (19) can be written as 
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Note that from Eqs. (20) and (21), and T^ ^ can be interpreted as the strain and 

kinetic energy respectively of elements associated with design variable x^. 

Thus, given p and y , the rate of change of eigenvalues can be computed from 
the energies associated with this design variable. Since most general-purpose 
finite-element codes provide element strain energy as an output option, one way to 
implement (22) is to calculate ^ and ^ by summing strain energy and kinetic 

energy for all elements controlled by design variable x . 


In the above formulation, we have made use of the form of the stiffness and mass 
matrices, Eqs. (11), (12). Not all structural elements can fit into these models but 
some important cases do. Some of these are tabulated in Table 1. 


Using Eqs. (15) to (18), it is possible to derive explicit equations for 
eigenvector derivatives as well as second-order derivatives of eigenvalues in terms 
of energies associated with various design variables. These are quite tedious and 
have not been accomplished so far. In the following, we will discuss the special 
case of B. = y . . 

a i 


EIGENSOLUTION SENSITIVITIES FOR THE CASE p. = y. 

r i fl i 

we can use chain rules to rewrite sensitivity derivatives 
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( 26 ) 


Using Eqs. (26), (27) and replacing all x_, x_ in Eqs* (2) to (8) by a and a , and 

r s r o 

making use of the orthogonality properties of normal modes 

<\>l Mf. = 0 if £ * j 

as well as the linearity assumptions (Eqs. (15) and (16)), we can derive, after 
considerable algebraic manipulation, the following results: 
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It should be noted that can be considered as the "cross mode" strain 

energy, since it is the work done by the elastic force in jth mode (i.e. moving 

through displacement in the £th mode. Similarly, can be considered as the 

"cross mode" kinetic energy. Thus, eigensolution sensitivity derivatives can be 
computed readily when these energy terms become available. 
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SENSITIVITY DERIVATIVES FOR MODIFIED SYSTEMS 


In iterative analysis, we frequently require the eigensolution derivatives of a 
system different from the nominal design. In these situations, assumed mode 
reanalysis [6-7] appears to be very efficient. Let AK and AM denote the change to 
stiffness and mass matrices, respectively. Then, in term of a^, we have 


ND 

AK = I (of. - 1)K. 
i=l 1 


and 
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Following the development in Ref. 7, the eigensolution of the modified system 
can be computed approximately by solving the following reduced eigenvalue problem 


K q = X M q 


(38) 


where 
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where 4> is the truncated modal matrix of the original system, and 
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Once (37) is solved, the eigenvectors of the modified system (|i. , in terms of physical 
coordinates, can be completed from 

ill. = <t> q. (43) 

i i 

For modified systems, the eigenvalue derivatives, Eqs . (28) and (30) are still 
applicable except and are now defined by 

V *rj ' 5 (44) 
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and the eigenvector derivatives can be computed from 
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where 


aq £ n 

do T ” , Z C £rj q j 0 7 ) 

r j=l J J 

C, rj is as defined by (34) or (35) with . , T £rj defined by (44) and (45). 


DISCUSSION 


In Eq. (7), the second-order derivatives of eigenvalues are shown to be 
dependent on eigenvector derivatives. In the present formulation, we can compute 
2 

8 A^/8a r 8a g Using Eq. (34) without the need to compute eigenvector derivatives 

explicitly. Once the derivatives with respective to or’s are known, chain rules can 
be used to compute the derivatives with respective to design variables x's (Eqs. 
(23) to (25)). 


NUMERICAL EXAMPLE 


The assumed mode reanalysis sensitivity derivative formulation has been 
implemented in a program which post-processes MSC/NASTRAN generated data. The first- 
order sensitivity data have been applied to improve analytical model using measured 
modal data [8] as well as synthesis of structures with multiple frequency 
constraints [9]. Recently the second order derivatives of eigenvalues (Eq. (30)) 
has also been implemented. 

A ten-bar cantilever truss structure, Fig. 1, is used to test the program. The 
ten members are grouped into 4 design variables as indicated in Figure 1. Starting 

o 

with a uniform structure with cross sectional area 10 in for all design variables, 
the optimal design program described in Ref. 7 is used to mode the first two 
natural frequencies from 13.3 and 37.8 Hz to 16 and 39.3 Hz, respectively. This is 
accomplished by a sequential linear programming formulation [7,9]. At each 
intermediate design, the eigenvalue derivatives are computed using reanalysis 
formulations. Table 2 defines the design history. Specifically, the designs at 
iteration No. A-0 and B-0 are analyzed exactly using MSC/NASTRAN. Three iterations 
are shown after each exact analysis. The eigensolution at designs A-l to A-3 and 
B-l to B-3 are computed using assumed mode reanalysis formulations. Four modes are 
used in each case. The first two natural frequencies are tabulated in Table 3. 

Also shown in Table 3 are the corresponding exact frequencies. From Table 3, it can 
be seen that the accuracy in frequency of assumed mode reanalysis formulation is 
very good. Tables 4 to 7 summarize 3^ /3 , /3x 2 , 3^^ and 3^/3x 2 , 

respectively. The results of these tables indicate that the sensitivity derivatives 
of modified system can be predicted quite accurately using the assumed-mode 
reanalysis formulation. 
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CONCLUDING REMARKS 


General procedures for computing eigensolution sensitivity derivatives for a 
class of problems have been proposed in this paper. Detailed formulations have been 
carried out for a special case. It is shown that the eigenvalue derivative with a 
design variable can be computed from strain energy and kinetic energy for that design 
variable. Furthermore, when the cross mode energy terms are available, assumed mode 
method can be used for eigensolution as well as associated sensitivity reanalysis. 
This efficient formulation has proved to be very effective in synthesis of structures 
with multiple frequency constraints [7,9]. Additionally, the present approach can 
be implemented in a post-processor of any finite-element programs without the need 
to modify the source code. 

Since the current formulation provides an efficient approach for computing 
second-order eigenvalue derivatives, it would appear that a second-order method for 
structural optimization with frequency constraints could be implemented efficiently. 
Finally, in view of the success of the formulation for the special case of (3^ = 

further development for the general case of (3^ t seems to be warranted. 

SYMBOLS 


K 

M 


£rj 


r 

AK 

AM 

*£ 

N 


n 

ND 


system stiffness matrix 
system mass matrix 
eigenvector of the £th mode 

eigenvalue of the £th mode 

modal matrix of original system 
cross-mode strain energy 

cross-mode kinetic energy 

rth design variable 

modification in stiffness matrix 
modification in mass matrix 

eigenvector of the th mode of the modified systems 

number of dof of the system 
number of modes computed, n < N 
number of design variables 
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TABLE I. - STIFFNESS AND MASS EXPONENTS FOR TABLE II. - DESIGN HISTORY OF 

SEVERAL COMMON STRUCTURAL ELEMENTS TEN BAR TRUSS 


Element 

Design variable 

p 

Y 

Truss 

Cross 

i 

1 

Membrane 

Thickness 

i 

1 

Plate bending 

Thickness * 

3 

1 

Beam bending 

Cross-sectional area 

2 

1 

Beam bending 

Section area moment 




of inertia* 

1 

0.5 


Iteration 

No. 

X 1 

2 

X 3 

x 4 

A-0 

10.0 

10.0 

10.0 

10.0 

A- 1 

12.76 

8.88 

5.0 

5.0 

A-2 

10.39 

8.15 

5.0 

5.0 

A-3 

10.44 

8.29 

5.0 

5.0 

B-0 

10.44 

8.29 

5.0 

5.0 

B-l 

7.92 

7.80 

3.44 

3.44 

B-2 

7.47 

7.18 

2.90 

2.90 

B-3 

7.19 

7.0 

2.68 

2.68 


Circular cross section 


TABLE III. - COMPARISONS OF NATURAL FREQUENCIES 


Case 

No. 

fj(Hz) 


Error 

(%) 

f 2 (Hz) 


Error 

(%) 

Approximate 

Exact 

Approximate 

Exact 

A- 1 

16.88 

16.60 

1.71 

40 . 16 

39.03 

2.88 

A-2 

15.96 

15.78 

1.14 

39.85 

38.99 

2.19 

A-3 

15.99 

15.82 

1.10 

40.00 

39.11 

1.27 

B-l 

15.90 

15.89 

0.11 

40.23 

39.87 

0.88 

B-2 

15.99 

15.98 

0.03 

39.98 

39.44 

1.37 

B-3 

15.99 

15.99 

0.00 

40.00 

39.33 

1.71 


TABLE IV. - COMPARISON OF a^/DXj TABLE V. - COMPARISON OF 3A 1 /3x 2 


Case 

No. 

3\/3x 

3A/3x 

Error 

a) 

Eq. (28) 

Exact 

A-l 

216.7 

294.6 

26.5 

A-2 

208.8 

269.9 

22.5 

A-3 

202.3 

262.7 

23.5 

B-l 

203.8 

205.2 

0.7 

B-2 

203.8 

222.9 

8.6 

B-3 

196.0 

221.0 

11.3 


Case 

No. 

3A/MX 

3A./MX 

Error 

(%) 

Eq. (28) 

Exact 

A-l 

418.9 

454.0 

7.7 

A-2 

494.2 

509.2 

2.9 

A-3 

495.8 

505.9 

1.9 

B-l 

557.2 

578.3 

3.7 

B-2 

607.3 

609.2 

0.3 

B-3 

635.8 

620.0 

2.6 
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Modulus of elasticity = 10^ ksi 

o 

Material density =0.1 lbm/in 

Design variable x^ = 1,2, 3, 4 

x = 5,6 

x^ = member 7,8,9,10 
Figure 1. Ten Member Cantilever Truss 



